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Recently S. Galam and A. Mauger [Phys. Rev. E 56, 
322 (1997)] proposed an approximant which relates the bond 
and the site percolation threshold for a particular lattice. 
Their formula is based on a fit to exact and simulation results 
obtained earlier for different periodic and aperiodic lattices. 
However, the numerical result for an aperiodic dodecagonal 
lattice does not agree well with the proposed formula. I 
present here new and more precise data for this and other 
aperiodic lattices. The previously published value for the do- 
decagonal lattice is confirmed. The reason for the deviation 
from the Galam and Mauger approximant is discussed. 
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In a recent article |^] S. Galam and A. Mauger pro- 
posed a formula which relates the site and bond perco- 
lation thresholds for a given lattice. This formula is a 
modification of previous one [|| designed to predict the 
percolation thresholds for a variety of periodic lattices. 
In the new Galam and Mauger approximant (NGMA) the 
(mean) coordination number of a lattice (z) is replaced 
by effective one {zeff)- 

Pc=PO [{d-l){Zeff-l]-''d'' 

where po, a and b are free parameters and the (countable) 
set of values for them divide lattices in separate classes. 
The formula provides a connection between the site and 
the bond percolation threshold in a particular class. So, 
if (say) the bond percolation threshold is known for a 
lattice (and one guesses somehow to which class belongs 
this lattice) then one can determine z^ff. After replacing 
Po, a and b with their values for site percolation (within 
the same class) one could obtain a value for the site per- 
colation threshold. In the same manner one could start 
from a site percolation threshold and to get the bond 
threshold. 

The NGMA was tested on simulation and exact values 
for 17 lattices and among them six lattices based on ape- 
riodic tilings. The values for two variants of the octag- 
onal and dodecagonal aperiodic lattices were extracted 
from my computer simulation results It appeared 
that some of these values do not fit well to the NGMA. 
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The authors of |i| supposed that the value for bond per- 
colation on the ferro variant of the dodecagonal lattice 
should be corrected. Instead of my result pc = 0.495 
they suggest that the correct value should be 0.475 - 
so it would fit much better to the theoretical curve (the 
NGMA). Since the bond percolation results given in my 
paper were stated to be preliminary, it is reasonable to 
be suspicious about them. 

Here I present new and more precise numerical data 
for these thresholds as well as for the bond percolation 
thresholds for the other aperiodic lattices used for testing 
the NGMA - see the left part of Table |. 

These values were obtained by computer simulation 
analogous to that in The lattice sizes were much 
larger here: up to size 500 x 500 bond lengths (the 
longer bond) and the computational efforts were be- 
tween 100 and 1000 times larger for the different lat- 
tices. For this work I used other pseudo random num- 
ber generators (PRNG). For the most part of simulation 
the drand48() PRNG from the standard GnuC distri- 
bution was used. I used also the ran2 generator from 
Q and the recently discussed |^ four tap shift-register 
generator R(471, 1586, 6988, 9689). The results coincided 
within less than 1/10 of a percent in the mean of the 
cumulative Gaussian distribution which was used for a 
fitting function. The fits just with a polynomial of 3-th 
degree gave almost the same results for the places (the 
values of p) where the spanning probability curve reaches 
a value of 1/2. One could pay attention also, that the 
value for the non-modified Penrose lattice conforms the 
earlier results. 

So I could claim (within the rigor of a Monte Carlo 
simulation) that the bond percolation threshold for the 
"ferro" modification of the dodecagonal aperiodic lattice 
is indeed 0.495 within error bars of ±0.001. In fact a 
more precise value could be extracted from the data on 
Fig.l pc = 0.4950 ± 0.0005. The error bars would be 
even smaller if one disregard the shift in the spanning 
probability curve for L = 300. 

Now one could draw the conclusion that deviations of 
this percolation threshold from the NGMA do exist and 
one has either to accept that the approximant is good 
within a larger interval of deviations: \pc — Pc\ = A < 
0.014 or to search a way to improve the approximant 
itself. 

The reason for the discrepancy could be suggested in 
the line of thinking which S. van der Marck has presented 
recently B for site percolation. He paid attention that 
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for certain lattices there are sites positioned on specific 
positions which do not contribute to percolation at all — 
e.g. in 2D these are the "star sites" inside a triangle. 

Indeed for bond percolation such entirely irrelevant 
bonds are impossible. However, one can show that some 
bonds contribute much less to the short range connectiv- 
ity than the others. This is the case of the questioned 
aperiodic dodecagonal lattice with "ferro" bonds. This 
lattice is a modification of the original dodecagonal lat- 
tice constructed by Socolar, which corresponds to an ape- 
riodic tiling of plane with hexagon, square and a rhom- 
bus. The "ferro" modifications consists in adding new 
bonds along the short diagonals of the rhombuses. Now 
one can show that these new bonds contribute less to the 
short range percolative connectivity. 

An estimation for this contribution could be done by 
stating a percolation problem for "one rhombus system" . 
We have five bonds arranged in a rhombus where the 
fifth bond is along the short diagonal. The question is 
what is the probability for a bond to be a "cutting" one 
for percolation between sites on the acute angles (i.e. 
what is the probability that the extraction of this bond 
destroys the connection between the two sites). One can 
easily see that this probability is 2p^{l — p)^ (the two 
zig-zag paths) for the ferro bond, which gives « 1/16 for 
p 1/2 (close to the questioned threshold value). The 
same value (near to p — 1/2) is « 3/16 for the other 
bonds [4p3(l -p)2 +p'^{l -p)^ +p^{l -p)]. 

Of course, things change when one considers percola- 
tion between the sites connected directly with the "ferro" 
bond. The cutting probability (again near p = 1/2) for 
this bond is already « 9/32 [p{l — p^)^] and the cut- 
ting probability of each one of the other bonds is « 3/32 

[p^(l-p)2(l+p)]. 

It appears that the chance to be a cutting bond is the 
same for the two type of bonds if both pairs of opposite 
vertices of the rhombus are considered . But one has to 
mention (e.g. on Fig. Ig in [Q): there are only one or 
two outgoing bonds from the obtuse angle vertices, while 
the acute vortices have 2 to 4 outgoing bonds. 

So if one considers the "spreading of connectivity" one 
more step beyond the two diagonals of the rhombus, the 
contribution of the "ferro" bond is less significant (by a 
factor of two approximately) . 

In contrast to the above, the same analysis for central- 
force rigidity percolation on the same lattices [^| shows 
that the "ferro" bond is always cutting bond within the 
questioned rhombus. There is a simple approximant for 
the rigidity thresholds - so called Maxwell approximant 
[0,1). It is just twice the other mean- field like approxi- 
mant: the Scher and Zallen prediction for (connectivity) 
bond percolation thresholds in two dimensions {pc ~ 2/2;, 
where z is the mean coordination number). One can see 
from the Table || that the estimates for rigidity percola- 
tion models, where each bond has (almost) equal signifi- 
cance for spreading the rigidity, agree with the Maxwell 
approximant extremely well. 
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FIG. 1. The frequency of occurrence of spanning cluster 
S(p,L) for dodecagonal lattices with different size(L). Each 
data point is averaged over 1000 realizations. 



TABLE 1. New numerical results for the bond percolation 
thresholds of questioned aperiodic lattices. The values in right 
are the rigidity (bond) percolation thresholds (pJJ) taken from 
[8]. Some of the lattices are always floppy (a.f.) in view of 
the rigidity percolation model. One can see that the lattices 
problematic for the NGMA have pc larger than the half of the 
respective rigidity threshold. The error estimates are put in 
brackets and concern the underlined digit. 
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Penrose 


0.4767(5) 
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Penrose (f) 


0.429(1) 


0.420... 


0.836(2) 


0.840... 


Octagonal 


0.478(3) 


1/2 


a.f. 
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Octagonal(f) 


0.402(5) 


0.387... 


0.769(2) 


0.774... 


Dodecagonal 


0.538(1) 


0.551... 


a.f. 


1.102... 


Dodecagonal(f) 


0.4950(5) 


0.468... 


0.938(1) 


0.937... 
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